Isospin equilibration in multi-fragmentation processes and dynamical correlations 
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The asymptotic time derivative of the total dipole signal is proposed as an useful observable 
to study Isospin equilibration phenomenon in multi-fragmentation processes. The study proceeds 
through the investigation of the 40 CZ + 28 Si system at 40 MeV/nucleon by means of semiclassical 
microscopic many-body calculations based on the CoMD-II model. In particular, the study has been 
developed to describe charge/mass equilibration processes involving the gas and liquid "phases" of 
the total system formed during the early stage of a collision. Through the investigation of dynamical 
many-body correlations, it is also shown how the proposed observable is rather sensitive to different 
parameterizations of the isospin dependent interaction. 



I. INTRODUCTION 

An interesting subject related to Heavy Ions Isospin 
physics p. J is the process leading to the equilibration 
of the charge/mass ratio between the main partners of 
the reaction as well described in Ref.[2]. The so called 
"isospin diffusion" phenomenon is the relevant mech- 
anism acting between the reaction partners in binary 
processes HI 03 E]. In particular, in the collision of 
the 124 and 112 Tin isotopes at 50 MeV/nucleon [3], 
evidence of partial equilibrium in the charge/mass ratios 
of the quasi-projectile and quasi-target has been deduced 
through the study of the iso-scaling parameters related to 
the isotopic distributions. In this case dynamical calcula- 
tions based on the Boltzmann-Uehling-Uhlenbeck model 
[7] show that the degree of equilibration depends on the 
behavior of the symmetry potential U T as a function of 
the density. The analysis of the experimental data in 
this kind of studies is based on the linear relation be- 
tween the iso-scaling parameter and the relative neutron 
excess (3 of the emitting sources (typical in several sta- 
tistical models). It is also assumed that both quantities 
weakly depend on the secondary statistical decay pro- 
cesses. In general, as discussed in ref. [5], this last condi- 
tion can be affected by the fragment excitation energies, 
(or temperatures) and by the distinctive features of the 
models used to simulate the secondary decay stage. In 
this work we want to extend the study of the isospin 
equilibration processes, looking at the whole system, by 
using the following quantity V*(i) = J2fl°i The sum 
on the index i is performed on all the Z tot protons of the 
system. V*(t) corresponds, apart from the elementary 
charge e, to the time derivative of the total dipole of the 
system. The velocities vf are computed in the center of 
mass (cm.) reference frame. Several studies were based 
on this dynamical variable to describe pre-equilibrium 7- 
ray emission (see Refs. [9j [TOj [11] and references therein 
). Various reasons suggest us to use the same variable to 
also describe isospin equilibration processes. 

- (i) After the pre-equilibrium stage, starting from the 
time t pre , when a second stage characterized by an av- 
erage isotropic emission of the secondary sources (sta- 
tistical equilibrium) takes place, the ensemble average 



of V (t) satisfies the following relation: v(t pre ) = V(t> 

t pre ) = \^ 1 10] - The average value of this dynamical vari- 
able at t pre is in fact invariant with respect to statistical 

processes and therefore the value of V is determined only 
by the complex dynamics which characterizes the early 
stage of the collision, when fast changes of the average 

nuclear density are expected. In particular, I? can be 
expressed as a function of the charge Z, mass A, aver- 
age multiplicity friz, a and the mean momentum (P)z,A 
of the detected particles having charge Z and mass A in 
the generic event: 
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C -U is the correlation function between the multiplic- 
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ity and the mean momentum. This correlation function 
plays a key role for the invariance property and therefore 
requires for an event by event analysis in which many- 
body correlations can not be neglected. 

For symmetry reasons, lies on the reaction plane. 
It is directly linked with a weighted mean of the 
charge/mass ratio, as Eq.(l) suggests. It also takes into 
account the average isospin flow direction through the 

momenta (P) z A . 

-(ii) In the general case , we find attractive the follow- 
ing decomposition: = \?q + \ + V*gl where 1? G 

and V l are the average dipolar signals associated to the 
gas "phase" (light charged particles) and to the "liquid" 
part [TJ], corresponding to the motion of the produced 

heavy fragments. The signal \^gl is instead associated 
to the relative motion of the two "phases" . By supposing, 
for simplicity, that the gas " phase" is formed by neutrons 

and protons, can be further decomposed as: 

V = - — V r H V cm ,LG + Vr,L 

(3) 
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In the above expression the first term represents the con- 
tribution related to the neutron-proton relative motion 
of the gas "phase" expressed through the relative ve- 
locity lf^ p ; the second term concerns the relative ve- 
locity ~T? C m,LG between the centers of mass of the "liq- 
uid" complex and the "gas"; the last term represents 
the contribution produced by the relative motion of the 
fragments. A similar expression can be obtained includ- 
ing others light particles in the gas "phase" as , for ex- 
ample, the Intermediate Mass Fragments (IMF). From 
this decomposition we can see how the isospin equilibra- 
tion condition (V? = 0), for the total system, requires a 
very delicate balance which depends on the average neu- 
tron excess (/3 = T ) of the produced "liquid drops" 
/3l, on the one associated to the gas "phase" fto, and 
on the relative velocities between the different parts. 
To enlighten the role played by some of the terms re- 
ported in Eq.(3), we can discuss the idealized decay of 
a charge/mass asymmetric source through neutron and 
proton emission (or the case in which the liquid drops are 

produced through a statistical mechanism V V,l = 0). 
Moreover, for simplicity, we can consider uncorrelated 
fluctuations between the velocities, masses and neutron 
excesses. In absence of pre-equilibrium emission or for 
identical colliding nuclei the second term of Eq.(3) is 
zero (lfcTn,LG = 0), and the isospin equilibration requires 
a neutron or proton gas "phase" or absence of relative 
neutron-proton motion. For non- identical colliding nu- 
clei, if pre-equilibrium emission exists, then Jf C m,LG 7^ 0- 
In this case, if (3q ^ Pli due, for example, to the isospin 
"distillation" phenomenon, the first term has to be nec- 
essarily different from zero and it will contribute to the 
neutron-proton differential flow (see also Sect. III). There- 
fore, according to our description, the understanding of 
the isospin equilibration process for the total system re- 
quires the gas "phase" contribution to be taken into ac- 
count. This term can be regarded as a kind of "dissi- 
pation" with respect to the system formed by the liq- 
uid part. In this work, as an example, we will discuss 
the results obtained through the Constrained Molecu- 
lar Dynamics-II approach (CoMD-II) [13j[14] applied to 
the charge/mass asymmetric system 40 C7 + 28 Si at 40 
MeV/nucleon. The study is performed by using different 
options for the symmetry potential term . 

Before to show the results of our calculations, in the 
following section we briefly recall the way in which the 
isospin interaction is introduced in CoMD-II model. 



II. SYMMETRY INTERACTION AND 
CORRELATIONS 

According to the results obtained in Ref. [15], in the 
simple case of large and compact systems, the isospin 
dependent part of the interaction is expressed, in the so 



called Non-Local (N.L.), approximation as: 

UkL - = ^ A2F ' [{l + \ a - a ' )f32 -\ a] (4) 
ld 2 a 

a = 4 W 1/3=0 (5) 

/3 4 terms are neglected in the previous expression, a = 
a{p) represents the correlation coefficient related to the 
difference in the dynamics of the n-p couples with re- 
spect to the n-n and p-p ones. It is evaluated at [3 — 
and describes the main effect of the many-body corre- 
lations on the iso-vectorial interaction, a depends on 
the average overlap integral per couples of nucleons p 
which reflects the degree of compression. The coeffi- 
cient a sym = 72MeV determines the strength of the iso- 
vectorial interaction at the ground state density ( g.s. ). 
F' is a form factor which modulates the changes of the 
iso-vectorial interaction as a function of the average over- 
lap integral S. In particular for the Stiff 1 option we use 
F' = Scj ] S +s , for the Stiff2 case F' = 1 and for the Soft 

option F' — (^c^-) 1 / 2 . From eq.(4) we note that the 
isospin forces if treated in a self-consistent many-body 
approach generates, beyond the (3 2 dependent potential, 
also another iso-vectorial density dependent term, not 
proportional to /3 2 , but proportional to the degree of cor- 
relation a. As discussed in [13 [TB], at small asymmetries, 
this term determines the high sensitivity of the experi- 
mental observable to the different functional form of F' . 
The limiting case corresponding to vanishing values for 
the correlation a represents the so called Iso-vectorial 
Mean Field Approximation (I.M.F.A.). In this case the 
average overlap integrals per couple of nucleons related 
to neutron-neutron, proton-proton and neutron-proton 
interactions have the same values and the iso-vectorial 
interactions generate only the usual symmetry potential 
term which is proportional to /3 2 (see eq. 4). A compar- 
ison between results obtained by full CoMD-II calcula- 
tions and the one derived in the I.M.F.A. limit will be 
discussed in the next section. 



III. CALCULATION RESULTS 

Now we discuss the results concerning the isospin 
equilibration process for the 40 C7 + 28 Si system at 40 
MeV/nucleon. For this collision, as an example, in Fig. 
1 we show the average total dipolar signals evaluated 
through CoMD-II calculations along the z beam direc- 
tion V and along the impact parameter direction x, 
V , respectively. The reference frame is the cm. one. 
The impact parameter b is equal to 3 fm, in panels (a) 
and (c), and 1.5 fm in panels (b) and (d). In Fig. 1(a) 
and Fig. 1(c), the average dipolar signals are shown for 
the first 150 fm/c. Different lines refer to different iso- 
vectorial potentials. The isospin independent compress- 
ibility is equal to 220 MeV, according to Ref. [13]. In the 
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FIG. 1: Average dipolar signals V for b=3 fm along the z 
direction are plotted as a function of time in the intervals 0- 
150 fm/c (panel (a)) and 100-400 fm/c (panel (b)). Different 
lines indicate different options for the isospin interaction (see 
the text). Panels (c) and (d) display the average dipolar sig- 
nals along the x direction for b=1.5 fm and in the same time 
intervals like panels (a) and (b) respectively. 
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FIG. 2: Asymptotic values V and V evaluated for different 
value of the b r parameter (see the text). Different symbols 
indicate different options for the iso-vectorial interaction. 



first 150 fm/ c we can see that in all the cases wide oscilla- 
tions exist. They are responsible for the pre-equilibrium 
7-rays emission [5J [TU]. The damped oscillations con- 
verge towards smaller and constant values (within the 
uncertainty of the statistics associated to the ensemble 
average procedure). This can be seen in Fig. 1(b) and 
Fig. 1(d) in which the dynamical evolution is followed 
from 100 fm/c up to 300 fm/c. The time interval in 
which the stationary behavior is reached is related to 
the life time of the coherent dipolar collective mode and 
it is strictly linked to the average time for the forma- 
tion of the main fragments and pre-equilibrium emission. 
For the asymptotic components, global results are shown 
in Fig. 2. They concern different impact parameters and 
different interaction options. In particular, in Fig. 3(a), 
we show for b=3 fm the charge distribution evaluated 
after 650 fm/c for the Stiff2 option. The shape of the 
distribution clearly show that at this impact parameter 
the mechanisms evolve essentially throughout a multi- 
fragmentation process. Under the same conditions, we 



show in Fig. 3(b) and in Fig. 3(c) the ratio Rql = jytj 

between z asymptotic components of the dipolar signal 
related to the light particles (first 2 terms of eq.(3)) and 
the one associated to the two biggest fragment (last term 
of eq. (3) ) . Rq l is plotted as a function of the reduced im- 
pact parameter b r = (b max ~ 7.5fm). In Fig. 3(b) we 
can see that for central and mid-peripheral collisions the 
"gas" component (pre-equilibrium contribution) is com- 
parable or larger than the one associated to the "liquid" 
one. Fig. 3(c) shows that the relevance of the "gas" light 
particle signal rapidly decreases with the increasing of 
the impact parameter. 

As discussed in the introductory section when the col- 
lision partners are not identical nuclei the isospin equi- 
libration process produce a neutron-proton differential 
flow contribution if the neutron and proton " gases" have 
different cm. velocities. These pre-conditions are veri- 
fied for the studied collision. For b = 3 fm and for the 
Stiff 1 and Soft options in Fig. 3(d) we show the neutron- 
proton differential flow F np as a function of the parti- 
cles rapidity, y, normalized to the projectile one ybeam- 
The rapidity values are evaluated in the cm. reference 
frame. From the figure we can see that, on average, the 
neutron-proton transversal velocity has a negative value. 
This reflects the "bending" of the relative motion be- 
tween the cm. of the emitted neutrons and protons, 
through the half-plane opposite to the impact parameter 
direction. These results can be compared with the cal- 
culations displayed in Fig. 3(e) obtained by subtracting, 
event by event, the cm. relative neutron-proton motion 
related to the "gas" phase . As can be seen, similarly 
to the case of identical nuclei, this correction restores 
(within the errors associated to the statistics of simula- 
tions) the almost specular of the F np behavior with re- 
spect the rapidity axes. The correction acts also along 
the beam direction. The final results also show an en- 
hanced sensitivity to the options related to the isospin 
potential. 

In the following we want to discuss in some detail the 
sensitivity of the dipolar signal to different options con- 
cerning the iso-vectorial interaction and different approx- 
imation scheme. In Fig. 4(a) we show as a function of the 
reduced impact parameter b r = (b max ~ 7.5fm) the 

asymptotic values of V* and V 2 . The arrow indicates the 
direction of increasing impact parameters corresponding 
to the marked points. Different symbols represent dif- 
ferent options. In the region of the bending of the lines, 
around 6 = 3 fm, we observe the greater sensitivity to 
the different options. This result is particularly evident 
by studying the ratios R = V /V . In Fig. 4(b) we in 
fact show the relative change r = between couples of 
different options. We can see that for b r less than about 
0.6 large changes are predicted according to the differ- 
ent form factor shapes. This impact parameter region is 
clearly dominated by large overlap between projectile and 
target nuclei which gives rise to processes changing from 
incomplete fusion reactions to IMF production. The re- 
gion of intermediate impact parameters show the higher 
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FIG. 3: (a) charge Z distribution of the reaction products ob- 
tained at 650 fm/c for the Stiff2 option and 6=3 fm. In panels 

I v z I 

(b) and (c), the Rgl = j^fj ratio (see the text) is plotted as 

a function of the reduced impact parameter. Panel (d) dis- 
plays neutron-proton differential flow F np as a function of the 
cm. rapidity. Different symbols refers to different options de- 
scribing the iso-vectorial interaction, (e) the same quantities 
reported in the panel (d) are plotted after the correction for 
the relative velocity associated to the cm. of the neutron an 
protons " gases" . 
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FIG. 4: -(a) V is plotted as a function of the corresponding 
V value for different b r values and different options (differ- 
ent symbols). The arrow indicates the direction of increasing 
impact parameters, -(b) relative changes r for the ratio R 
(see the text) evaluated for different couples of options as a 
function of the b r reduced impact parameter, -(c) Values of r 
evaluated for the Stiff l-StifF2 options as a function of b r are 
plotted in the case of full CoMD-II calculations and in the 
case of the I.M.F.A. approximation (see the text). The lines 
which join the points are meant only to guide the eye through 
the shown trend. 



sensitivity when the mechanism evolves with respect to 
the essentially binary processes which take place at the 
higher impact parameters. For b r greater than 0.6, in 
fact, the sensitivity is strongly reduced. 

In particular, according to what previously observed 
(see for example point (ii)), we have evaluated the par- 
tial contributions V L and V L related to the two main 
fragments. As an example, for 6 = 3 fm and for the 
Stiff2 option, the "liquid" asymptotic values are V L — 
—0.120c and V L = 0.162c while the total contributions 
are V = 0.044c and V = —0.027c. It results therefore 
that the contributions carried by the two main fragments 
only partially contribute to the isospin equilibration pro- 
cess. The remanent part ("gas"), which in this case we 
have associated to particles and to the IMF, generates 
a term with opposite sign and similar strength for both 
directions. It contributes in a decisive way to the global 
equilibration process. In particular, as an example, for 
b = 3 fm, in Fig 4(b) we show with the star symbol, the 
sensitivity parameter r evaluated by changing the option 
from Stiffl to Stiff2 and by only taking into account the 
two main fragments contributions. As we can see, the 
partial contribution shows a rather reduced sensitivity 
to the different options as compared to the case obtained 
by using the complete information of the system. 

Finally, in the following, we show the role of the correla- 
tion coefficient a, introduced in Sec. II, into determine the 
sensitivity of the investigated observable to the density 
dependence of the iso-vectorial interaction. For this aim, 
in Fig. 4(c), we compare, for different impact parameters, 
the values of r obtained for the Stiff2-Stiffl options with 
the ones obtained in the I.M.F.A. limit. Fig. 4(c) clearly 
shows that in the I.M.F.A. case, at the investigated en- 
ergies, the sensitivity of our investigated phenomenon to 
the behavior of the symmetry interaction is rather re- 
duced. The I.M.F.A. also strongly affects the values of 

T^. In particular, independently from the used options, 
it produces for 6,.=0.4, values of \V | about four times 
larger than the ones obtained with full COMD-II calcula- 
tions indicating a reduced capacity to obtain the isospin 
equilibration along the z direction. 



IV. CONCLUSIVE REMARKS 

In summary, in this work the isospin equilibration pro- 
cess for the asymmetric charge/mass system 40 CI + 28 Si 
at 40 MeV/nucleon has been investigated by studying 
the ensemble average of the time derivative of the to- 
tal dipole evaluated through CoMD-II calculations. 
Some general properties of this quantity have been dis- 
cussed. In particular, it allows to generalize the definition 
of isospin equilibration also in complex reactions evolving 
through multi-fragmentation processes. CoMD-II calcu- 
lations show that the asymptotic values of \? for these 
processes are quite sensitive to different symmetry po- 



5 



tential options; moreover, in central and mid-peripheral, 
the dipolar contribution associated to the pre-equilibrium 
emission of charged particles is relevant to determine the 

value of V and the related sensitivity to different density 



dependent form factors. CoMD-II calculations performed 
in the so called I.M.F.A. scheme also highlights the fun- 
damental role played by the many-body correlations in 
the study of the isospin equilibration processes. 
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